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Abstract. Supernova neutrino flavor transitions during the shock wave propagation 
are known to encode relevant information not only about the matter density profile 
but also about unknown neutrino properties, such as the mass hierarchy (normal 
or inverted) and the mixing angle ^13. While previous studies have focussed on 
"deterministic" density profiles, we investigate the effect of possible stochastic matter 
density fluctuations in the wake of supernova shock waves. In particular, we 
study the impact of small-scale fluctuations on the electron (anti)neutrino survival 
probability, and on the observable spectra of inverse-beta-decay events in future 
water-Cherenkov detectors. We find that such fiuctuations, even with relatively small 
amplitudes, can have significant damping effects on the fiavor transition pattern, and 
can partly erase the shock-wave imprint on the observable time spectra, especially for 
sin^ 6*13 > 0(10-3). 



PACS numbers: 14.60.Pq, 97.60.Bw, 02.50.Ey 
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1. Introduction 

Future observations of supernova (SN) neutrinos in underground detectors represent a 
subject of intensive investigation in astroparticle pliysics. In this context, matter effects 
associated to neutrino flavor transitions in tfie SN envelope have been widely studied as 
a unique tool to probe, at the same time, neutrino properties and supernova astrophysics 
(see, e.g. [1-4] for recent reviews). In particular, the effects of supernova shock waves 
on neutrino flavor transitions in matter [5] are gaining increasing attention in the recent 
literature [6-11]. Indeed, for a few seconds after the SN core bounce, the strong time 
dependence of the shock-wave profile can leave peculiar signatures on the time structure 
of the neutrino events, which could be monitored in large, real-time detectors [12-14]. 

So far, studies of supernova neutrino flavor transitions during the shock-wave 
propagation have been based on "deterministic" matter density profiles, assumed to 
be known (or at least knowable, in principle) both in time and in radial dependence. 
However, stochastic density fluctuations and inhomogeneities, of various magnitudes and 
correlation lengths, may reasonably arise in the wake of a shock front. Possible causes of 
these inhomogeneities include microscopic fluctuations in the nascent neutron star [15] 
and large-scale fluctuations between the proto-neutron star and the supernova envelope 
due to hydrodynamical instabilities [16-18]. At early times (<1 s after bounce), post- 
shock convection overturns can also produce large density anisotropics. A supernova 
neutrino "beam" traveling to the Earth might thus experience stochastic matter effects 
while traversing the stellar envelope. 

Concerning neutrino oscillations, the phenomenology of possible stochastic matter 
density fluctuations has been investigated in several contexts, with emphasis on general 
properties [19-22] and on the solution to the solar neutrino problem [23-28]. Supernova 
neutrinos have specifically been considered in relatively few cases [29], and only for a 
static profile. In general, it is found that the typical effect of random fluctuations on 
neutrino oscillations is to wash out the phase information (if any) and to damp the flavor 
transition pattern in the energy domain. It seems thus worthwhile to revisit this topic 
in the context of dynamic SN density profiles, which offer a complementary handle in 
the time domain. J Moreover, stochastic density fluctuations are expected to arise more 
easily in the wake of a SN shock wave than in a relatively static stellar environment. 

In this work we try to explore quantitatively the entangled effects of shock waves 
and of possible stochastic fluctuations (behind the shock front) on supernova neutrino 
flavor transition probabilities and observables. Unfortunately, current core-collapse SN 
simulations do not yet offer a clear input for such phenomena, both because the details 
of the explosion mechanism are not well understood yet, and because current computer 
resources do not allow to resolve density variations at scales smaller than, say, 0(10) 
km. Therefore, some assumptions and simpliflcations are unavoidable. In particular, we 
shall limit ourselves to small-scale and small-amplitude fluctuations, which are deflnitely 
not excluded by current simulations, and which allow a simple perturbative approach. 

I This possibility was mentioned in passing in the seminal Ref. [5]. 
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Our work is organized as follows. In Section 2 we introduce the notation 
for the neutrino mass and mixing parameters and for the electron (anti) neutrino 
survival probability Pee characterizing supernova neutrino transitions. In Section 3 we 
parametrize the "fluctuating" shock-wave profile and discuss an approximate analytic 
expression for the survival probability Pee, which neatly includes the damping effect 
induced by density fluctuations. In particular, the analytical derivation of Pee from 
the neutrino master equation for the density matrix, together with a comparison with 
representative numerical solutions, are given in the Appendix. We find that, for 
sin^ ^13 ^O(10~^), small-scale stochastic fiuctuations can suppress the imprint of the 
shock wave on the flavor transition pattern in the time domain. In Section 4 we 
discuss a speciflc experimental application for the case of antineutrinos in inverted 
hierarchy, by calculating the effects of fluctuations on positron event spectra observable 
in water- Cherenkov detectors through inverse beta decay. Conclusions and prospects 
for further developments are given in Section 5. The bottom line of our work is that 
some "smearing" of shock-wave signatures in neutrino observables is to be expected, if 
stochastic matter fluctuations are present in the wake of the shock front. Improvements 
in supernova hydrodynamical simulations will greatly help to clearify this issue. 



2. Three-neutrino mixing framework 

In this section we set the notation for three-neutrino mixing and for the survival 
probability of electron (anti)neutrinos Pee- In the current "standard" Su scenario, the 
dominant parameters relevant to flavor transitions in supernovae are the largest squared 
mass difference Am^ and the mixing angles 9i2 and ^13 (see, e.g., [30,31]). In numerical 
examples, we take as reference values (close to their best flts [32]) 

Am^ ~ 2.4 X 10~^ eV^ , (1) 
sin^ 612 ~ 0.3 . (2) 

The sign of Am^ distinguishes the cases of normal hierarchy (NH: +Am^) and inverted 
hierarchy (IH: — Am^). Concerning the mixing parameter sin^^i3, we shall use 
representative values below the current upper limits (sin^^i3 < few x 10~^ [33]). 

The smallest squared mass difference (5m^ ~ 8 x 10~^ eV^ [32]) is such that 
(5m^/Am^ <C 1; together with the smallness of sin^^i3, this fact guarantees, to a very 
good approximation, the factorization of the Su dynamics into a "low" (L) and a "high" 
(H) 2v subsystem. In other words, the electron (anti) neutrino survival probability Pge, 
up to small terms of 0(sin^ ^^13, 5m^/Am^), can be expressed as 

Pee — Pee ' Pee i (3) 

where P^^ and P^^ represent effective electron (anti) neutrino survival probabilities in 
the 2z/ subsystems (see [30,31,34-36] and references therein). 

Neglecting Earth matter effects (not included in this work), one has [8,30,31] 

pL^\ sin^6'i2 (forz/, any hierarchy), 
I cos^6'i2 (for any hierarchy). 
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and 

pH ^ f P+ (for in NH or V in IH), 

ee-'ypH ^fQj, p -^^ j^fj i/ in IH), ^ ^ 

where are defined below. 

In general, P^ depends on ^13, on the hierarchy, and on the wavenumber kn in the 
H subsystem, 

±kH = ±Am^/2E (+ for NH, - for IH), (6) 
and on the (anti)neutrino potential in matter [37], 

±V{x) = ±V2Gf Ne{x) (+ for u, - for V) (7) 

where is the electron density at the supernova radius x. 

Strong matter effects arc generally expected when the "level crossing condition" 

±kH ^ ±V{xc) (for u in NH or V in IH) , (8) 

is satisfied at (more than) one point Xc (see [31]). The crossing condition requires equal 
signs for V and k^, so it is not realized for V in NH or i/ in IH. The physics of matter 
effects can be encoded in terms of the so-called crossing probability Pc [30, 38] 

p ^ f Pcikn, sin' ^13, V) (for u in NH or V in IH), 
" \-0 (forI7inNHori/inIH), 

where Pc 7^ (=0) defines the case of nonadiabatic (adiabatic) matter transitions. § 

In the standard case (with no stochastic fluctuations) the survival probability P^ 
at the Earth (averaged over many oscillation cycles) is related to the crossing probability 
Pc through the well-known Parke's formula [44], 

Pf ~ i + (^i - Pc) cos 2^13 cos 2ef^{xo) (no fluctuations) , (10) 
where Xo is the neutrino production point, and 

cos2^~f3(x) = cos2e,,TV{x)/kn ^^^^ 

'(cos 2^13 T V{x)/kH f + (sin 2^13)2 



defines the mixing angle in matter. In Eq. (11), the upper and lower signs correspond 
to the upper and lower signs of P^ in Eq. (5). Due to the high matter density at the 
origin {\V{xQ)/kH\ ^ 1), in Eq. (10) one can simply put 

cos2^±3(xo) ~Tl . (12) 

A final remark is in order. In the presence of density fiuctuations with small 

amplitude, one docs not expect that the crossing probability Pc is significantly 
perturbed, being related to a "local" nonadiabatic effect. The relation between P^ and 

§ Concerning the early literature on two-level crossings, it is amusing to note that the classic Landau- 
Zener-Stiickelberg results for [39-41] were similarly obtained by Majorana [42] in the context of 
spin-flip transitions in variable magnetic fields, see the discussion in [43]. Curiously, this contribution 
by Majorana (well known in atomic physics) is largely ignored in the neutrino literature. 
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Pc in Eq. (10) is instead expected to change, being related to the "global" propagation 
history within the supernova matter. This intuitive picture, which enters in the 
analytical generalization of Eq. (10) discussed in the next section, is confirmed by 
numerical calculations (see the Appendix). Similarly, one does not expect that small- 
amplitude fluctuations can ruin the effective L®H factorization, which relies only on the 
smallness of ^13 and of bvr? j tyrr? . This expectation is also supported by the smallness 
of fluctuation effects within the L subsystem in our framework (see the last paragraph 
in Sec. 3.3). A full numerical conflrmation of the L®H factorization in the presence of 
stochastic effects, however, would involve rather heavy calculations of exact 3i/ solutions, 
which are beyond the scope of this work. 



3. Effects of stochcistic matter fluctuations on SN neutrino transitions 

In this Section we parametrize the SN stochastic density fluctuations through some 
simplifying assumptions (Sec. 3.1), characterize their effects through an analytical 
presciption for the calculation of the probabihty (Sec. 3.2), and discuss the time 
dependence of P^ in the presence of fluctuations with increasing amplitude (Sec. 3.3). 

3.1. Shock-wave density profile and stochastic fluctuations 

We make the reasonable assumption that stochastic fluctuations arise only afler the 
passage of the shock-wave. Behind the shock front, fluctuations can be described as 
fractional (random) variations ^(x) of an "average" neutrino potential Vq, 

V{x) = V^{x) + 5V{x) = V^{x)[l + i{x)] , (13) 

with vanishing mean {i{x)) = 0, and with variance = {^{x)^). In the absence of 
further information, we simply assume that ^ is both constant and small [0(few%)].|| 

We also assume that the fluctuations arise at relatively small length scales Lq, as 
compared with the neutrino oscillation length in matter A^. In the region where matter 
effects are relevant (F ~ /cj^), this condition reads: 

Lo < A„ ~ - — — — 
Kh sm 26^13 

, / E \ (2x 10-3 cV^X / 0.2 \ , 

- ( i(nfev) ( A,„^ j ■ (1*) 

In all numerical examples we shall assume a representative value 

Lo ~ 10 km . (15) 

Notice that this scale is of the same order as the radial hydrodynamical fluctuations 
described in [15], as well as of the ncTitrinosphere radius [1], which sets the transverse 
size of the SN neutrino beam in the Earth direction. 

II In a diflFcrcnt context (Sun density profile), the combination of all available solar and KaniLAND 
data sets an upper limit ^ < 5% at 70 % CL for the solar density fluctuations, assuming a fluctuation 
length scale Lq — 10 km [28]. 
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u potential V (+ fluctuations with ^ = 4 % ) 




Figure 1. Simplified radial profiles of the neutrino potential V{x) at different post- 
bounce times t (1, 2, 4, 8, and 16 s). The thicker curves behind the shock front 
represent fluctuations of amplitude ^ < 4% of the local matter potential V. Upper 
panel: Forward shock only. Lower panel: Forward plus reverse shock. In both cases, 
the static (pre-bounce) profile (i < s) is also shown. The band within dashed lines 
marks the region where SN matter effects are potentially important {V ~ fc//, with 
E = 2-60 MeV). 



The small-scale assumption in Eqs. (14) and (15) implies that the autocorrelation 
function of the random field ^{x) can be written as an effective delta-function [20-23] 

{axinx2)) = 2LoeS{x,-X2) . (16) 

In a nutshell, fluctuations between points much closer (much farther) than one neutrino 
oscillation length in matter are assumed to be fully correlated (totally uncorrelated) . 
Given Eq. (16), the master equation for the neutrino flavor evolution is determined (see 
the Appendix). 

Figure 1 shows the (fluctuating) shock-wave profiles used in this work. We 
take the unperturbated shock- wave potentials Vo{x) from our previous semplified 
parametrizations in [8, 12]. In the region behind the shock front, the thicker curves 
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represent a band of fluctuations with fractional deviation ^ < 4%. The upper panel 
refers to the neutrino potential V in the presence of the forward shock only. In the 
lower panel we have taken into account also a possible reverse shock, characterized by a 
smaller density jump at the front [9]. In both panels, we also show the band spanned by 
the neutrino wavenumber = Am'^/2E for E e [2, 60] MeV. Notice that the (crossing) 
condition for large matter effects {V ~ kn) is basically unperturbed by fluctuations as 
small as those in Fig. 1. 

3.2. Analytical results 

As anticipated in Sec. 2, the survival probability in the H subsytem is modified by 
matter density fluctuations. In general, the calculation of requires that the stochastic 
master equation (see Appendix) is solved numerically. Nevertheless, as shown in [24] and 
described in detail in the Appendix, small-scale and small- amplitude fluctuations allow 
to use perturbativc tcchniqiics which, already at first order, provide rather accTirate 
analytical approximations. In the context of SN neutrinos, the analytical prescription 
reads 

P±-l'^{l- P) e"^^ cos 2^13 , (17) 

which generahzes [24] Parke's formula (10) [in the conditions of Eq. (12)] through a 
simple exponential damping factor e"'"^'^, whose exponent is given by 

r±= r V{x) sm^ 2ef^{x)dx , (18) 

where 

= eVo^(x)Lo , (19) 

and the mixing angle in matter {Of-^) is defined as in Eq. (11). The domain of the 
damping integral V± is the region behind the forward shock front position x^, where 
stochastic fiuctuations are assumed to arise. For strong damping {T± ^ 1) one gets the 
limit P^ ^1/2, corresponding to a sort of complete "flavor depolarization," where the 
two effective v states in the H subsystem arc democratically mixed. 

All stochastic matter effects are embedded in the damping factor e~'"±, while the 
crossing probability Pc = Pc(^_ff, sin^ ^^is, V (a;)) in Eq. (17) is basically unperturbed.^ 
In particular, P^ can be evaluated analytically by considering the ordered product 
of crossing matrices along the average (nonmonotonic) supernova shock-wave density 
proflle Vo{x), as previously discussed in [8,12]. Therefore, apart from the numerical 
evaluation of the integral in Eq. (18), P^ can be computed through convenient analytical 
approximations. 

^ This fact is also verified a posteriori through the very good agreement of the analytical recipe in 

Eq. (17) with the numerical solution, as shown in the Appendix. The check is notrivial, since the 
accuracy of the analytical approximation in Eq. (17), discussed in [24] in the context of the (static 
and monotonic) solar matter profile, cannot be taken for granted in the (dynamic and nonmonotic) 
supernova profile. 
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3.3. Analysis of the survival probability in the H subsystem, and comments on 
in the L subsystem 

In this Section we describe the behaviour of as a function of time, in the presence 
of stochastic fluctuations with increasing amphtude ^. The case of no fluctuations is 
recovered for ^ = 0. 

Figure 2 shows the variations of the survival probability P^ (relevant for v in NH 
or V in IH) , 

P+-\-i(^- P) cos 2^13 , (20) 

in the time interval t G [1, 13] s, for ^ = (solid curves), ^ = 2% (dashed curves) and 
^ = 4% (dotted curves). The neutrino energy is fixed at £^ = 30 MeV, while sin^^ia 
ranges from 10~^ (upper panels) to 10~^ (lower panels). We include both the case of 
forward shock only (left panels) and of forward plus reverse shock (right panels). 

In the absence of fluctuations, the function P^{t) in Fig. 2 shows strong, 
nonmonotonic variations: A clear signatTirc of the strongly nonadiabatic flavor 
transitions along the shock front, operative in the time window t G~ [3, 9] s (as discussed 
at length in [8,12]). In the presence of stochastic fluctuations behind the shock front, 
however, the variations are partly suppressed, as a result of the "flavor depolarization" 
effect. As the amplitude ^ increases, the survival probability P^ gets closer to the 
saturation value 1/2. The effect of fluctuations is relevant in the same time window 
where shock effects are operative, since the occurrence of the condition y ~ /cjj (in 
multiple points Xc) leads to sin^ 29^^^ ~ 1 and thus to a "large" damping factor [see 
Eq. (18)] around Xc- At earlier or later times, the enhancement of ^13 in matter is more 
modest and so are fluctuation effects. 

In Fig. 2, the smearing effect of fluctuations appears to be particularly dangerous for 
the identification of shock-wave signatures: For instance, for ^ = 4%, the double-peak 
structure associated to the forward+reversc shock (right panels) may be washed out, 
becoming more similar to the case of forward shock only (left panels). A "confusion 
scenario" might arise, limiting the potential to monitor the shock-wave propagation 
through neutrino flavor transitions. 

Figure 2 also shows that, in general, damping effects are suppressed at small values 
of ^13 (say, sin^^is^lO-^). This fact can be understood by expanding the damping 
integral [Eq. (18)] around the crossing points x^., where the integrand is locally maximal. 
At first order in x — Xc, and using the fact that V ku around Xc, one gets then 

(21) 

which shows that the damping factor is basically proportional to ^13. In conclusion, 
for values of sin^ ^13 <; O(10~^), small-scale stochastic fluctuations with a fractional 
amplitude of a few percent might significantly suppress shock-wave effects on the electron 
neutrino survival probability P_?. For smaller values of sin^ 6*13, fluctuation effects appear 



r+ ~ TT sin 2^13 ^P( 



d\nV 
dx 
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forward shock forword + reverse shock 




2 4 6 8 10 12 2 4 6 8 10 12 14 

time (s) time (s) 

Figure 2. Probability function -Pf^(i) (relevant for v in NH or V in IH) a.t E — 30 MeV, 
for four representative values of sin^ ^13, in the presence of density fluctuations with 
fractional amplitude ^ equal to 0% (solid curves), 2% (dashed curves) and 4% (dotted 
curves). Left panels: forward shock only. Right panels: forward plus reverse shock. 



to be less dangerous, but unfortunately the overall effect of shock waves on P_|^ is also 
small. 

Figure 3 shows the time evolution of the survival probability (relevant for z/ in 
IH or V in NH, where Pc = 0), 

P^~^(l + e-r- cos 2^13) , (22) 
in the same format of Fig. 2. In the absence of matter fluctuations (^ = 0), it is r_ = 



Damping of supernova neutrino transitions in stochastic shock-wave density profiles 10 



and ~ 1 trivially (being cos 2^13 ~ 1 in all cases shown in Fig. 3) . In the presence of 
fluctuations, the damping term e~^- tends to lower (down to the plateau value 1/2 
for large T ). However, the relative variations of the function in Fig. 3 are not as 
strong as for P^ in Fig. 2, since the integrand in r_ is never "resonant", as it happens 
instead for r+. For small 6*13, it is simply 

r_ ^ sin^ 29,s r , (23) 

Jxo (1 + v{x)/kHy 

namely, the damping factor r_ decreases with ~ ^13) and its effect vanishes more rapidly 
than for r_|_. This fact explains why, in Fig. 3, fluctuation effects are visible only at larger 
values of ^13 (upper panels, sin^ 613 — 10^^) as compared with Fig. 2. The absence of 
level crossing effects in P^ strongly reduces the sensitivity to the details of the shock- 
wave profile, so that for ^ > there is hardly any difference between the left and right 
upper panels in Fig. 3. 

From the previous discussion, stochastic density fluctuations appear to be most 
relevant in the cases where P^ (rather than P^) is involved, i.e., neutrinos in normal 
hierarchy or antineutrinos in inverted hierarchy. For this reason, in the next section we 
shall examine the impact of these results in a phenomenologically relevant case {V in 
IH) , by calculating observable spectra of positrons induced by u^. through inverse beta 
decay. 

Finally, we comment about the effects of fluctuations on the survival probability 
in the "low" (L) subsystem. In the presence of density fluctuations, this probability 
can be calculated as 

P,^,^iTie-r±(^) cos 2^12 , (24) 

which is analogous to Eq. (17) for the adiabatic case Pc = (the adiabaticity of the L 
subsystem is discussed in [8]). The damping integral in the above equations reads 

T±{L)= rV{x) sin^ 2e^2{x)dx , (25) 

Jxo 

where the deflnition of 6'^2 is analogous to that in Eq. (11), but in terms of 612 and of 
ki — Sm'^/2E. Since ki is a factor ~ 30 smaller than kn, one has tipically V/kL ^ I 
for X < Xs, and thus sin^ 29i2{x) <^ 1 in Eq. (25), except possibly at late times {t ^ 10 s, 
see Fig. 1), where the crossing condition V ^ k^ can be realized, leading to sin^^i2 — 1 
for neutrinos. However, at late times the neutrino luminosity is also small. Therefore, 
the damping effects of T±{L) are either small or suppressed by a low luminosity, and 
can be safely neglected for the purposes of our work. One can simply take e~^^^^^ ~ 1 
in Eq. (24), thus recovering Eq. (4). 

4. Observable positron spectra from SN V in inverted hierarchy 

As an application of the previous results, we study the effects of stochastic matter 
fluctuations on the observable positron time spectra detectable through the inverse 
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Figure 3. As in Fig. 2, but for the probability function P^{t) (relevant for ly in IH 
or F in NH) . 



beta-decay reaction 

'T'e + P ^ n + e'^ , (26) 

which represents the main SN neutrino detection channel in present [45] and planned 
Cherenkov detectors [46-49], as well as in scintillation detectors [50-52]. Therefore, 
in this Section we focus on electron antineutrinos, in the relevant case of inverted 
hierarchy,^ where shock-wave (plus fluctuation) effects become potentially observable 

An analogous study (not presented here) could be performed for the neutrino channel in normal 
hierarchy, in the context of z^e-sensitive detectors, such as liquid argon time proportional chambers [53]. 
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through the H subsystem dynamics, 

Pee-COS^^lsPf . (27) 

Positron event rates are calculated for a reference 0.4 Mton water- Cherenkov 
detector, with neutrino spectra at the supernova source, interaction cross sections, and 
detection parameters fixed as in our previous work [12], to which the reader is referred 
for further details. Here we just mention that in [12] two representative "low" and 
"high" positron energy bins (at E'pos = 20 ± 5 and 45 ± 5 MeV, respectively) were 
identified as good relative tracers of shock-wave effects. Indeed, as discussed in [12] 
(for no fluctuations), the time spectra at relatively high energy (-Epos = 45 ± 5 MeV) 
show strong signatures of the shock- wave propagation, especially for increasing sin^ ^13, 
through nonmonotonic time variations of the event rates. Conversely, the shock-induced 
time variations in the bin Ep^s = 20 ± 5 MeV are significantly smaller; at such energy, in 
fact, the initial electron and non-electron antineutrino fiuxes (as used in our calculation 
input) happen to be approximatively equal, and the effects of fiavor transitions largely 
cancel. 

Figure 4 shows absolute time spectra of events within the two reference positron 
energy bins (-Epos = 20 ± 5 MeV and -Epos = 45 ± 5 MeV) for four representative 
values of sin^^^i3, and for both ^ = (no fluctuations, solid histograms) and ^ = 4% 
(dashed histograms). The left and right panels refer to the case of forward shock only 
and of forward plus reverse shock, respectively. In the presence of stochastic density 
fluctuations (^ = 4% in Fig. 4), the spectra for the low-energy positron bin are basically 
unaffected, since all flavor transition effects (modifled or not by fluctuations) are small. 
The time spectra for the high-energy positron bin (-Epos = 45 ± 5 MeV) are instead 
signiflcantly smoothed out by fluctuation effects for sin^ ^13 ^0(10"^), as expected from 
the discussion of Fig. 2. Therefore, even at small amplitude = 4%), fluctuations 
can suppress the shock imprint on the time spectra, and make them qualitatively 
similar to those in normal hierarchy (where they are expected a priori to be smooth 
and monotonic). 

As suggested in [12] (for the case of no fluctuations), the relative flavor transition 
effects at low and high energy are best displayed by showing the ratio of the spectra 
at 20 ± 5 MeV and 45 ± 5 MeV, where the flrst spectrum acts as a "normalization" 
factor. Figure 5 shows such ratio (in inverted hierarchy) for the same values of ^ and 
sin^ ^13 of Fig. 4. In the absence of stochastic fluctuations (^ = 0), we recover the fact 
[12] that the spectral ratio can track the time variations of the electron antineutrino 
survival probability , as evident from a visual comparison of the solid curves (^ = 0) 
in Figs. 2 and 5.* Unfortunately, stochastic density fluctuations (dashed curves with 
{ = 4 % in Fig. 5) destroy this nice correspondence for sin^ 6*13 ^ O{10^^). The spectral 
ratio becomes smoother, and the "peaks and valleys" induced by the shock wave on P^ 
disappear. The damping effect implies a net loss of information about the shock wave 

* The analysis in [12] actually refers to the crossing probability Pc- However, in the absence of 
fluctuations, the smallness of ^13 implies that « Pc, see Eq. (10). 
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Figure 4. Absolute time spectra of positron events induced by in IH in a 0.4 Mton 
water-Cherenkov detector, in the presence of forward shock (left panels) and forward 
plus reverse shock (right panels) . Four representative values of sin^ ^13 are considered. 
The solid histograms refer to the case of no fluctuation (^ = 0), while the dashed ones 
refer to the case of fluctuations with ^ = 4%. In each panel, the upper (lower) couple of 
histograms refers to the positron energy bin i?pos = 20± 5 MeV (-Epos = 45± 5 MeV). 



in the observable neutrino spectra. 
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Figure 5. Time dependence of the ratio of events between the energy bins Epos = 
20 ± 5 MeV and ii^pos = 45 ± 5 MeV in the case of toward shock only (left panels) and 
forward plus reverse shock (right panels). The ratio refers to in inverted hierarchy 
for no fluctuation (^ = 0, solid curves) and a fluctuation amplitude ^ — 4% of the local 
matter potential (dashed curves). 



5. Summary and conclusions 

Stochastic density fluctuations in supernova matter may produce observable effects on 
the supernova neutrino signal. In this context, previous works have been focussed 
on static and monotonic proflles. In this paper, we have investigated the case of 
time-dependent and non-monotonic proflles, embedding forward (plus reverse) shock 
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propagation, as suggested by recent SN numerical simulations. In the hypothesis 
of small-scale (Lq ~ 0(10 km)) and small-amplitude < fcw%) fluctuations, we 
have discussed an analytical recipe to evaluate the SN electron (anti) neutrino survival 
probability Pgo which accounts for both standard matter transitions and additional 
"flavor-depolarization" effects induced by the fluctuations, and which accurately 
reproduces the numerical results (see the Appendix). 

We find that, stochastic fiuctuations — possibly arising after the shock front 
passage — may significantly suppress the imprint of SN shock waves on Pee in the time 
domain, the more the larger sin^ ^13. For decreasing sin^ ^13, fiuctuation effects decrease, 
but shock-wave signatures also become less pronounced. An application to observable 
time spectra of positron events in large water-Cherenkov detectors shows that, for the 
phenomenologically relevant case of inverted hierarchy, the time spectra can easily loose 
any imprint of the shock wave for sin^ ^13 ^ (9(10^^). Therefore, in the presence of 
stochastic fluctuations, it might be difficult to "monitor" the shock-wave in real time 
in a future galactic SN explosion, or to find unmistakable signatures of inverted mass 
hierarchy effects. 

In this context, future work might include a "fiuctuation analysis" of the SN 1987A 
neutrino events [54] and of the expected supernova relic neutrino spectrum [55-57]. 
For instance, some studies of the SN 1987A in the case of inverted hierarchy [58] 
found a tension between data and theory in the case of adiabatic transitions (i.e., 
sin^ ^^13 ^ 10^'^) and inverted hierarchy, using a static SN profile. The effect of both 
nonadiabatic crossings and of damping behind the shock front might put these results 
in a different perspective. 

In any case, further studies will greatly benefit from improvements of numerical 
SN simulations, so as to follow the shock wave evolution with much higher space-time 
resolution than currently possible. A better theoretical understanding of the possibility 
and property of stochastic density fiuctuations in the SN envelope might also help to 
remove the simplifying assumptions adopted in this work. 
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Appendix A. Stochastic neutrino master equation: analytical and 
numerical solutions 

For the sake of completeness, in this Appendix we explicitly derive Eq. (17) (based 
on previous literature on the subject) and compare it with our numerical results for 
representative SN shock-wave cases. To avoid an excessive "±" notation, we consider 
here only the case of neutrinos in normal hierarchy. The cases of antineutrinos and 
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of inverted hierarchy are obtained by the replacements V — > —V and kn — —kn, 
respectively. 

In the presence of stochastic matter fluctuations, one must consider the neutrino 
evolution within a statistical ensemble of density profiles. Averaging over the ensemble 
generally produces a loss of the coherence in the system. In this context, the appropriate 
formalism involves the neutrino density matrix p. 

For a 2i/ mixing problem (as in the H subsystem), the density matrix in the flavor 
basis {i/e, ^aY^ 1 where Va is any linear combination of and can be expressed in 
terms of a "polarization vector" P [59-62]: 

p=l(l + P-a), (A.l) 



where 



/ 2Re(z/>„,) \ 
21m(z/>,) 



(A.2) 



and a is the vector of Pauli matrices. The length |P| of the polarization vector 
(|Pp = Trp^) measures the degree of coherence of the system: |P| = 1 corresponds 
to a pure state, < |P| < 1 to a partially mixed state, and |P| = to a completely 
mixed state. 

With or without fluctuations, the Hamiltonian for neutrinos propagating in matter 
in the "high" subsystem is the usual one [37]: 

H=\-a^]^[kH sin 2^13] (Ti + ^ [V{x) - kn cos 2^13] (73 , (A.3) 

with 

h = (A;j^ sin 2^13, 0, \/(a;) - fc^^ cos 2^13) . (A.4) 

In the presence of stochastic matter density fluctuations, the potential has been 
expressed as V{x) = Vo{x)[l + ^{x}], with 5-correlated Gaussian fluctuations (see 
Sec. 3.1). In this case, the Liouville evolution of the fluctuations-averaged density matrix 
reads [20-22] ft 

= -i[Ho, p] - lp[(73, [as, p]] , (A.5) 

where Hq is the unperturbated Hamiltonian (for ^ = 0) and V is the same damping 
parameter as in Eq. (19). 

Equation (A.5) has the well-known Lindblad form for dissipative systems [63], where 
the flrst term in the right-hand side reproduces the standard Liouville equation for the 
density matrix, while the second term violates the conservation of Trp^ and allows 
transitions from pure to mixed states. In terms of the polarization vector, Eq. (A.5) 
assumes the neat form of a Bloch equation [60, 65] 

-^P = ho X P - VPt , (A.6) 
ax 

jj For neutrinos propagating in a stochastic matter, the density matrix [Eq. (A.l)] and the polarization 
vector [Eq. (A.2)] must be understood as averaged over the ensamble of density profiles. 
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where ho corresponds to Eq. (A. 4) with { = 0, and Pt = P — P3 is the "transverse" part 
of the polarization vector. In this picture, the first term in the right hand side represents 
the usual coherent "precession" of the polarization vector P around the vector ho [59- 
62,64], while the second (damping) term produces a loss of coherence, i.e., a shortening 
of the vector length |P| [60,65]. 

Using Eqs. (A.2)-(A.4), one can write Eq. (A. 6) in components [20-24] 



d_ 
dx 



( pA 

= K: , (A.7) 

[PsJ 

where 020 = cos 2^13 and S2e = sin 2^13. The (fluctuation-averaged) survival probability 
for electron neutrinos = (1, 0)-^ is then given by 

I + P3 



( pA 






-Vo{x) + kHC2e 





\ 


(pA 


P2 




Vq{x) - kHC2B 


-V{x) 


-kHS2e 




P2 


[PsJ 




[ 


kHS20 







[Ps) 



P^ = Tr[pz/V; 



(A.8) 



In general, the system of equations [Eq. (A.7)] has to be solved numerically, as done 
e.g. in [20,21,23]. However, in the case of small-amplitude fluctuations, the damping 
term V in Eq. (A.7) can be treated as a perturbation. In this way, it is possible to 
integrate analytically Eq. (A.7). Here we closely follow the derivation presented in [24]. 

We start by diagonalizing the /C matrix in Eq. (A.7) at 0-th order in P, estabihshing 
the "instantaneous diagonal basis" 



(A.9) 



(Pi] 




(pA 


P2 




P2 


[Psj 




[Ps) 



The matrix U — [ui, U2, Uq] diagonalizes /C as 
/Cd diag(Ai, A2, Ao) = UKU^ , 

where the three eigenvalues are 
Ao =0 , 



Ai,2= ±i|ho| = ±i^(A;ij cos 2^13- V)2 + (A;^^ sin 2^13)2 
and the corresponding eigenvectors are 

Uo 



Ul,2 



ho , 

1 

V2 



(ho X s 



(A.IO) 

(A.ll) 
(A.12) 

(A.13) 
(A.14) 



with ho = ho /I ho I, while s is an arbitrary unitary vector perpendicular to ho- By 
choosing s = (0, 1, 0), one obtains 



/ 



Uo = 



sin 2^J3 


— cos 29~{^ j 



COS 2^1+3 \ 



sin Idxz I 



I cos2^^3 \ 



sin 29 



(A.15) 



13 ) 
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while the matrix U takes the form 

1 „„„o^+ 1 



U 



^cos2e^3 ^cos2^^3 sin 2^1+3 
J_ _ J_ n 

^ sin 2^+3 ^ sin 2^+3 -cos 2^+ 



13 



where ^13 is the mixing angle in matter defined by 



sin 2^ 



sin 20 



13 



13 



cos 2^ 



^(cos2^i3 - V/kny + (sin 2^13)2 ' 
cos 2^13 - V/kn 



13 



^(cos2^i3-1^Aif)' + 

(sin 2^13)2 

In the diagonal basis, Eq. (A. 7) becomes 
ax ax 



(A.16) 

(A.17) 
(A.18) 

(A.19) 



We assume that the propagation is adiabatic {dUydx = 0), except near one 
"crossing point" Xc where V{xc) — kn (The generahzation to multiple crossings is 
straightforward). In this approximation, Eq. (A.19) has thus a formal solution: 

P(x) = U\x) S{xc, x) Q(xe) S{xo, Xc) U{xo) P(xo) , (A.20) 

where xq is the neutrino production point, and 



S{xi,X2) = diag 



(A.21) 



is the adiabatic evolution operator, while 0(xi,X2) = J^^ dy\ho{y)\, and the non- 
diagonal operator Q{xc) takes into account the nonadiabatic transition (level crossing) 
around the point Xc- 

By inserting the above expressions in Eq. (A. 8), and averaging over the oscillatory 
terms, one gets: 



_ 1 + [Q(^c)]33 COS 2g+3(x) COS 2^+3(^0) 

+ ~ 2 



(A.22) 



The comparison of Eq. (A.22) with Eq. (10) leads to the identification [Q(a;c)]33 = 
1 — 2Pc, where Pc is the crossing probability between the two mass eigenstates in matter. 
At 0-th order in T> one thus recovers the usual (no-fiuctuation) Parke's formula: 



1 



1 



2 , ^--Pe) COS 2^+3(0;) COS 2^+3(xo) . 



(A.23) 

The first-order corrections to the eigenvalues [Eqs. (A.11)-(A.12)] and eigenvectors 
[Eqs. (A.13)-(A.14)] are calculated through the standard perturbation theory as: 

AW = A„ + SXa = A„ - Vnigua , (A.24) 



Aa — Aft 



(A.25) 



where a, 6 = 0, 1, 2 and G — diag(l, 1, 0). Since u^^Ua = ^ab — (Ua)3(u6)3, one obtains 



A' 



13 



A;2 = ±i|ho| -v(i 



- sin=^ 26 
2 



13 



(A.26) 
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and 



f,{') _ r, , sin 2^+3 cos 2g+ ^ 

Uq — Uo + JJ TT— j S , 



^ (/) ^ sin 297n 



sin2^^3^ 

- cos 2^+3 Uo + ' U2,i 



(A.27) 



V2\ho 

By substituiting the first-order eigenvalues in Eq. (A. 20), the matrix S{xi,X2) can be 
rewritten as: 

i4>(xi,X2)+ f^^ dxSXi —i4>(xi,X2)+ dx S\2 +f^^dx5Xo 

e ■''"1 , e ■'^1 , e ■'^1 



S{xi,X2) = diag 



.(A.28) 



The expressions of 6\i^2 are irrelevant for our pourposes, since they disappear after 
averaging over the oscillating terms in P^. Notice that, in Eq. (A. 20), only the operator 
S is corrected at first order. The matrix ZY(xo) remains unchanged, since the high-density 
condition V{xo)/kH ^ 1 implies that sin26'i3(xo) ~ in Eq. (A.27). The matrix U{x) 
at the detection point is also unchanged, since ^13(2;) = ^13 and P = at the Earth. 

Finally, for small fluctuations, we assume that the operator Q{xc) (which only 
depends on the crossing condition V{xc) — kn) is unchanged. After averaging on the 
oscillating terms one then obtains 

e"^^ cos 2eUx) COS 2eUxo) , (A.29) 

where 

r+= rdyX',{y)= rdyV{y)sm^29Uy) . (A.30) 

Jxo Jxo 

In the case of detected supernova neutrinos, Oi^{x) corresponds to the vacuum 
value ^13, while the initial high-density condition implies cos2^i^(xo) — —1. With these 
positions, Eq. (A.29) gives the desired expression 

P^c^l-(l- p\ e-^+ cos 2^13 . (A.31) 



2 V2 

In order to check the reliability of the analytical approximation of P^ in Eq. (A.31), 
we have compared it with the results of a numerical (Runge-Kutta) evolution of 
neutrino master equation [Eq. (A. 7)] along representative fluctuating shock- wave 
profiles. Figure Al shows our calculation of P^{t) at fixed neutrino energy {E — 
30 MeV) for sin^^i3 = 10~^, and for two representative fiuctuation amplitudes, ^ — 2% 
(upper panels) and ^ = 4% (lower panels). We consider both the case of forward shock 
only (left panels) and forward plus reverse shock (right panels). It can be seen that 
the analytical calculations (solid curves) and the numerical Runge-Kutta calculations 
(dashed curves) are in very good agreement in the whole time interval t G [1, 13] s. 
The results of Figure Al (and of other checks that we have performed for different 
values of sin^ 6*13 and of the neutrino energy) show that the analytical calculation of P_[^ 
reproduces the numerical results with high accuracy. The same reassuring check has 
been performed also in the case of P^ (not shown). 
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Figure Al. Comparison of analytical and numerical calculations of the electron 
neutrino survival probability in the H subsystem (solid and dashed curves, 
respectively) for E = 30 MeV, sin^ ^13 = 10^^, and fluctuations with amplitude ^ = 2% 
(upper panels) and ^ = 4% (lower panels). The left panels refer to the case of forward 
shock only, while the right panels to the case of forward plus reverse shock. 



Damping of supernova neutrino transitions in stochastic shock-wave density profiles 21 
References 

[1] K. Kotake, K. Sato and K. Takahashi, "Explosion Mechanism, Neutrino Burst, and Gravitational 
Wave in Core-Collapse Supernovae," Kept. Prog. Phys. 69, 971 (2006) [astro-ph/0509456]. 

[2] G. G. RafFclt, "Supernova neutrino oscillations," in the Proceedings of Nobel Symposium 2004, 
Neutrino Physics (Haga Slott, Enkoping, Sweden, 2004), Phys. Scripta T121, 102 (2005) [hep- 
ph/0501049]. 

[3] A. Dighe, "Supernova neutrinos: Production, propagation and oscillations," in the Proceedings 

of Neutrino 2004, 21st International Conference on Neutrino Physics and Astrophysics (Paris, 

Prance, 14-19 Jun 2004), Nucl. Phys. Proc. Suppl. 143, 449 (2005) [hep-ph/0409268] 
[4] F. Cavanna, M. L. Costantini, O. Palamara and F. Vissani, "Neutrinos as astrophysical probes," 

Surveys High Encrg. Phys. 19, 35 (2004) [astro-ph/0311256]. 
[5] R. C. Schirato and G. M. Fuller, "Connection between supernova shocks, flavor transformation, 

and the neutrino signal," astro-ph/0205390 (unpublished). 
[6] K. Takahashi, K. Sato, H. E. Dalhed and J. R. Wilson, "Shock propagation and neutrino oscillation 

in supernova," Astropart. Phys. 20, 189 (2003) [astro-ph/0212195]. 
[7] C. Lunardini and A. Yu. Smirnov, "Probing the neutrino mass hierarchy and the 13-mixing with 

supernovae," JCAP 0306, 009 (2003) [hep-ph/0302033]. 
[8] G. L. Fogli, E. Lisi, A. Mirizzi, and D. Montanino, "Analysis of energy- and timc-dcpcndcnce of 

supernova shock effects on neutrino crossing probabilities," Phys. Rev. D 68, 033005 (2003) 

[hep-ph/0304056]. 

[9] R. Tomas, M. Kachelriefi, G. Raffelt, A. Dighe, H. T. Janka and L. Schcck, "Neutrino signatures of 
supernova shock and reverse shock propagation," JCAP 0409, 015 (2004) [astro-ph/0407132]. 

[10] S. Kawagoe, "Shock Wave Propagation in Prompt Supernova Explosion and the MSW Effect of 
Neutrino," talk at TAUP 2005, 9th International Conference on Topics in Astroparticle and 
Underground Physics (Zaragoza, Spain, 2005), available at www.unizar.es/taup2005/talks.htm. 

[11] B. Dasgupta and A. Dighe, "Phase effects in neutrino conversions during a supernova shock wave," 
hep-ph/0510219. 

[12] G. L. Fogli, E. Lisi, A. Mirizzi and D. Montanino, "Probing supernova shock waves and neutrino 
flavor transitions in next-generation water-Cherenkov detectors," JCAP 0504, 002 (2005) [hep- 
ph/0412046]. 

[13] V. Barger, P. Huber and D. Marfatia, "Supernova neutrinos can tell us the neutrino mass hierarchy 
independently of flux models," Phys. Lett. B 617, 167 (2005) [hcp-ph/0501184]. 

[14] S. Choubey, N. P. Harries and G. G. Ross, "Probing neutrino oscillations from supernovae shock 
waves via the IceCube detector," hep-ph/0605255. 

[15] W. Keil, H. T. Janka and E. MuUer, "Ledoux-convection in protoneutron stars: A clue to supernova 
nucleosynthesis?," Astrophys. Journ. 473, Llll (1996) [astro-ph/9610203]. 

[16] K. Kifonidis, T. Plewa, H. T. Janka and E. Mueller, "Non-spherical core collapse supernovae. I: 
Neutrino-driven convection, Rayleigh- Taylor instabilities, and the formation and propagation of 
metal clumps," Astron. Astrophys. 408, 621 (2003) [astro-ph/0302239]. 

[17] R. Buras, H. T. Janka, M. Rampp and K. Kifonidis, "Two-Dimensional Hydrodynamic Core- 
Collapse Supernova Simulations with Spectral Neutrino Transport II. Models for Different 
Progenitor Stars," astro-ph/0512189. 

[18] L. Scheck, K. Kifonidis, H. T. Janka and E. Mueller, "Multidimensional Supernova Simulations 
with Approximative Neutrino Transport I. Neutron Star Kicks and the Anisotropy of Neutrino- 
Driven Explosions in Two Spatial Dimensions," astro-ph/0601302. 

[19] R. F. Sawyer, "Neutrino oscillations in inhomogeneous matter," Phys. Rev. D 42, 3908 (1990). 

[20] F. N. Loreti and A. B. Balantekin, "Neutrino oscillations in noisy media," Phys. Rev. D 50, 4762 
(1994) [nucl-th/9406003]. 

[21] A. B. Balantekin, J. M. Fetter and F. N. Loreti, "The MSW effect in a fluctuating matter density," 
Phys. Rev. D 54, 3941 (1996) [astro-ph/9604061]. 



Damping of supernova neutrino transitions in stochastic shock-wave density profiles 22 

[22] F. Benatti and R. Florcanini, "Dissipativc neutrino oscillations in randomly fluctuating matter," 

Phys. Rev. D 71, 013003 (2005) [hep-ph/0412311]. 
[23] H. Nunokawa, A. Rossi, V. B. Semikoz and J. W. F. Valle, "The effect of random matter density 

perturbations on the MSW solution to the solar neutrino problem," Nucl. Phys. B 472, 495 

(1996) [hcp-ph/9602307]. 

[24] C. P. Burgess and D. Michaud, "Neutrino propagation in a fluctuating sun," Annals Phys. 256, 1 

(1997) [hep-ph/9606295]. 

[25] A. A. Bykov, M. C. Gonzalez-Garcia, C. Pena-Garay, V. Y. Popov and V. B. Semikoz, "MSW 
solutions to the solar neutrino problem in presence of noisy matter density fluctuations," hep- 
ph/0005244. 

[26] C. Burgess, N. S. Dzhalilov, M. Maltoni, T. I. Rashba, V. B. Semikoz, M. Tortola and J. W. F. VaUe, 
"Large mixing angle oscillations as a probe of the deep solar interior," Astrophys. J. 588, L65 
(2003) [hep-ph/0209094]. 

[27] M. M. Guzzo, P. C. de Holanda and N. Reggiani, "Large mixing angle solution to the solar 
neutrino problem and random matter density perturbations," Phys. Lett. B 569, 45 (2003) 

[hep-ph/0303203]. 

[28] A. B. Balantekin and H. Yuksel, "Do the KamLAND and solar neutrino data rule out solar density 

fluctuations?," Phys. Rev. D 68, 013006 (2003) [hep-ph/0303169]. 
[29] F. N. Loreti, Y. Z. Qian, G. M. Fuller and A. B. Balantekin, "Effects of random density fluctuations 

on matter enhanced neutrino flavor transitions in supernovae and implications for supernova 

dynamics and nucleosynthesis," Phys. Rev. D 52, 6664 (1995) [astro-ph/9508106]. 
[30] A. S. Dighe and A. Yu. Smirnov, "Identifying the neutrino mass spectrum from the neutrino burst 

from a supernova," Phys. Rev. D 62, 033007 (2000) [hcp-ph/9907423]. 
[31] G. L. Fogli, E. Lisi, D. Montanino and A. Palazzo, "Supernova neutrino oscillations: A simple 

analytical approach," Phys. Rev. D 65, 073008 (2002) [Erratum-ibid. D 66, 039901 (2002)] 

[hep-ph/0111199]. 

[32] G. L. Fogli, E. Lisi, A. Marrone and A. Palazzo, "Global analysis of three-flavor neutrino masses 
and mixings," hep-ph/0506083; to appear in Prog. Part. Nucl. Phys (2006). 

[33] M. ApoUonio et al, "Search for neutrino oscillations on a long base-line at the CHOOZ nuclear 
power station," Eur. Phys. .J. C 27, 331 (2003) [hcp-ex/0301017]. 

[34] T. K. Kuo and J. T. Pantaleone, "Supernova Neutrinos And Their Oscillations," Phys. Rev. D 
37, 298 (1988). 

[35] S. P. Mikhcyev and A. Yu. Smirnov, "Resonant Neutrino Oscillations In Matter," Prog. Part. 

Nucl. Phys. 23, 41 (1989). 
[36] T. K. Kuo and J. T. Pantaleone, "Neutrino Oscillations In Matter," Rev. Mod. Phys. 61, 937 
(1989). 

[37] L. Wolfenstcin, "Neutrino Oscillations In Matter," Phys. Rev. D 17, 2369 (1978); S. P. Mikhcev 
and A. Yu. Smirnov, "Resonance Enhancement Of Oscillations In Matter And Solar Neutrino 
Spectroscopy," Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 913 (1985)]. 

[38] S. T. Petcov, "Exact Analytic Description Of Two Neutrino Oscillations In Matter With 
Exponentially Varying Density," Phys. Lett. B 200, 373 (1988). 

[39] L. D. Landau, "Zur Theorie der Energieiibertragung bei Stossen" ("On the theory of the energy 
transfer in collisions"), Phys. Z. Sowjetunion 1, 88 (1932); "Zur Theorie der Energieiibertragung, 
II" ("On the theory of the energy transfer, 11"), Phys. Z. Sowjetunion 2, 46 (1932). 

[40] C. Zener, "Nonadiabatic Crossing Of Energy Levels," Proc. Roy. Soc. Lond. A 137, 696 (1932). 

[41] E. C. Stiickelberg, "Theorie der unelastichen Stosse zwischen Atomen" ("Theory of the inelastic 
collisions between atoms"), Helv. Phys. Acta 5, 369 (1932). 

[42] E. Majorana, "Atomi orientati in campo magnetico variabile" ("Oriented atoms in variable 
magnetic field"), Nuovo Cimento 9, 43 (1932). 

[43] F. di Giacomo and E. E. Nikitin, " Majorana formula and the Landau-Zener-Stiickelberg treatment 
of the avoided crossing problem," Ups. Fiz. Nauk 175 (5), 545 (2005) [Physics-Uspekhi 48 (5), 



Damping of supernova neutrino transitions in stochastic shock-wave density profiles 23 



515 (2005)]. 

[44] S. J. Parke, "Nonadiabatic Level Crossing in Resonant Neutrino Oscillations," Phys. Rev. Lett. 
57, 1275 (1986). 

Y. Fukuda et al, "The Super-Kamiokande detector," Nucl. Instrum. Meth. A 501, 418 (2003). 

UNO official homepage: ale.physics.sunysb.edu 

C. K. Jung, "Feasibihty of a next generation underground water Cherenkov detector: UNO," 
hep-ex/0005046. 

K. Nakamura, "Hypcr-Kamiokande: A next generation water Cherenkov detector," Int. J. Mod. 

Phys. A 18, 4053 (2003). 
L. Mosca, "Status of the project of a Large International Underground Laboratory at 
Prejus," talk at the Villars CERN/SPSC Meeting (Villars, Switzerland, 2004), available at 
nuspp.in2p3.fr/Frcjus. See also J. E. Campagnc, M. Maltoni, M. Mezzetto and T. Schwetz, 
"Physics potential of the CERN-MEMPHYS neutrino oscillation project," hep-ph/0603172. 
L. Cadonati, F. P. Calaprice and M. C. Chen, "Supernova neutrino detection in Borexino," 

Astropart. Phys. 16, 361 (2002) [hep-ph/0012082]. 
M. Aglictta et ai, "The most powerful scintillator supernovae detector: LVD," Nuovo Cim. A 105, 
1793 (1992). 

L. Oberauer, "Low energy neutrino physics after SNO and KamLAND," Mod. Phys. Lett. 

A 19, 337 (2004) [hcp-ph/0402162]; L. Oberauer, "LENA," talk at NO-VE 2006, 3rd 
International Workshop on Neutrino Oscillations in Venice (Venice, Italy, 2006), available at 
http://axpd24.pd.in£n.it /NO-VE2006 /talks/NOVE_Oberauer.ppt. 
I. Gil-Botella and A. Rubbia, "Decoupling supernova and neutrino oscillation physics with LAr 

TPC detectors," JCAP 0408, 001 (2004) [hcp-ph/0404151]. 
A. Mirizzi and G. G. Raffelt, "New analysis of the SN 1987A neutrinos with a flexible spectral 

shape," Phys. Rev. D 72, 063001 (2005) [astro-ph/0508612]. 
M. Malek et al. [Super-Kamiokande Collaboration], "Search for supernova relic neutrinos at Super- 
Kamiokande," Phys. Rev. Lett. 90, 061101 (2003) [hcp-cx/0209028]. 

G. L. Fogli, E. Lisi, A. Mirizzi and D. Montanino, "Three-generation flavor transitions and decays 
of supernova relic neutrinos," Phys. Rev. D 70, 013001 (2004) [hep-ph/0401227]. 

H. Yukscl, S. Ando and J. F. Bcacom, "Direct measurement of supernova neutrino emission 
parameters with a gadolinium enhanced Super-Kamiokande detector," astro-ph/0509297. 

H. Minakata and H. Nunokawa, "Inverted hierarchy of neutrino masses disfavored by supernova 

1987A," Phys. Lett. B 504, 301 (2001) [hep-ph/0010240]. 
S. P. Mikhcycv and A. Yu. Smirnov, "Neutrino Oscillations In Matter With Varying Density," in 
the Proceedings of '86 Massive Neutrinos in Astrophysics and in Particle Physics, Sixth Moriond 
Workshop (Tignes, Prance, 1986), edited by O. Fackler and J. Tran Thanh Van (Editions 
Frontieres, Gif-sur-Yvcttc, 1986), p. 355. 
L. Stodolsky, "On The Treatment Of Neutrino Oscillations In A Thermal Environment," Phys. 
Rev. D 36, 2273 (1987). 

C. W. Kim, J. Kim and W. K. Szc, "On The Geometrical Representation Of Neutrino Oscillations 

In Vacuum And Matter," Phys. Rev. D 37, 1072 (1988). 
H. Duan, G. M. Fuller and Y. Z. Qian, "Collective neutrino flavor transformation in supernovae," 
astro-ph/0511275. 

G. Lindblad, "On The Generators Of Quantum Dynamical Semigroups," Commun. Math. Phys. 
48, 119 (1976). 

A. Yu. Smirnov, "Neutrino conversion and neutrino astrophysics," hcp-ph/9811296, in the 
Proceedings of Symposium on New Era in Neutrino Physics (Tokyo, Japan. 1998), published in 
"Tokyo 1998, New era in neutrino physics", edited by H. Minakata and O. Yasuda, Frontiers of 
Science Series n. 25 (Univ. Academy Press, Tokyo, Japan, 1999), p 1. 
[65] G. G. Raffelt, "Stars as Laboratories for Fundamental Physics," (Univ. of Chicago Press, Chicago, 
IL, 1996). 



